Phase diagrams of the 2D t — i! 



U Hubbard model from an extended mean field 
method 
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It is well-known from unrestricted Hartree-Fock computations that the 2D Hubbard model does 
not have homogeneous mean field states in significant regions of parameter space away from half 
filling. This is incompatible with standard mean field theory. We present a simple extension of 
the mean field method that avoids this problem. As in standard mean field theory, we restrict 
Hartree-Fock theory to simple translation invariant states describing antiferromagnetism (AF), fer- 
romagnetism (F) and paramagnetism (P), but we use an improved method to implement the doping 
constraint allowing us to detect when a phase separated state is energetically preferred, e.g. AF 
and F coexisting at the same time. We find that such mixed phases occur in significant parts of 
the phase diagrams, making them much richer than the ones from standard mean field theory. Our 
results for the 2D t — t' — U Hubbard model demonstrate the importance of band structure effects. 



PACS numbers: 71.10.Fd,05.70.Fh,75.50.Ee 



I. INTRODUCTION AND MAIN RESULTS 

Hubbard-type models in two dimensions have been fre- 
quently studied in the context of high temperature super- 
conductivity and other strongly correlated systemsi De- 
spite considerable efforts (for review see e.g. Ref.0) there 
is still need for simple methods that can contribute to the 
understanding of the complex behavior of such models. 
In this paper we study an extension of mean field (MF) 
theory which allows for the possibility of phase separated 
states, in addition to the usual MF states. We calculate 
full phase diagrams for the 2D t — t' — U Hubbard model 
which, to our knowledge, are not available in the litera- 
ture by other methods. 

MF theory offers several advantages compared to 
more complicated methods like unrestricted Hartree- 
Fock (HF) theory: It is easy to implement, not restricted 
to small system sizes, and can produce phase diagrams 
for Hubbard-type models with a limited computational 
effort. The disadvantage of standard MF theory is that it 
always predicts translation invariant states everywhere in 
the phase diagram, without giving any information about 
the stability with respect to fluctuations, or about the 
stability with respect to competing non-uniform states. 
In 2D Hubbard-type models these problems have severely 
restricted the usefulness of the MF approach*^ and the 
MF method is therefore not widely used. More correct 
methods indeed demonstrate that the qualitative features 
of the standard MF predictions are restricted to parts 
of the phase diagram, e.g., the antiferromagnetic (AF) 
phase at half filling. This suggests that the MF approach 
is unsatisfactory and motivates using more complicated 
methods. However, the more accurate theoretical meth- 
ods tend to be computationally demanding and therefore 
restricted to very small system sizes. 

In this paper we adopt and clarify the extended MF 



method in Refs.^laand use it to calculate phase diagrams 
of the 2D t — t' — U Hubbard model. This method is de- 
signed to overcome the limitation of only producing uni- 
form MF solutions, without increasing the computational 
effort. We use the standard mean field equationsj&il but 
we extend them by a method allowing us to detect pos- 
sible instabilities towards phase separation. 6 The phase 
diagrams we thus obtain are much richer than the ones 
obtained with conventional MF theory and no longer in 
contradiction with unrestricted HF results. In particu- 
lar, conventional MF theory for the 2D Hubbard model 
{if = 0) predicts an AF phase in a finite doping regime 
around half filling (see Fig. 3 in Ref. Q) which is known 
to be qualitatively wrong. The phase diagram from the 
extended MF theory is shown in Fig. 1. It shows that the 
AF phase exists only strictly at half filling, and at finite 
doping close to half filling no simple translation invariant 
state is thermodynamically stable, in agreement with un- 
restricted HF theory&2*i2iii While our method does not 
account for fluctuations or details of states which are not 
translational invariant, it allows to detect frustration in 
the sense of incompatibility between MF states and the 
doping constraint. Such frustration suggests interesting 
physical behavior to be explored by more sophisticated 
methods. Our theory should be useful also for other cases 
where no other methods are available. 

Our main results are the full phase diagrams for 2D 
t-t'-U Hubbard model for t' = and if = -0.35* 
in Figs. 1 and 2, respectively. They were obtained for 
a system size so large that they are practically identical 
with the thermodynamic limit. The phase diagrams are 
remarkably rich and very different from the correspond- 
ing results from standard MF theory: compare our Fig. 1 
with Fig. 3 in Ref. [3J and our Fig. 2 with Fig. I in Ref. 
4. Our results demonstrate that mixed phases are a typ- 
ical feature of 2D Hubbard-type models: as one changes 



zv 


• 

• 


i 


•-, 1 . • 

• • 




• 




• • 




• 




• • 




• 




• • 


15 


-• 
• 


F * 


• F *- 

• • 




• 




• • 




• 




• • 




• 




• • 


3 10 




* AF 

• • • 
• • • 


• m - 

• • • 
• • • 






• 

• mixed 


• 
mixed • 


5 




• 


• 






p 

• 


p 

• 






• 


• 


n 




• 


• 



-0.5 




p-1 



0.5 



20 



15 



3 10 



5 - 



mixed 



F 



* • • • 

• • • • 



mixed * 



-teW 



-0.5 0.5 

p-1 



FIG. 1: Phase diagram of the 2D Hubbard model as a func- 
tion of U and doping p for parameters t = 1 and t' — 0. 
We use Hartree-Fock theory restricted to ferromagnetic (F), 
antiferromagnetic (AF) and paramagnetic (P) states, and we 
find large mixed regimes where neither of these translational 
invariant states is thermodynamically stable. The results are 
for L = 60 and /3 = 1000 which is practically indistinguishable 
from the thermodynamic limit. 



doping one never goes directly from one MF phase to 
another, but there seems always a finite doping regime 
with a mixed phase in between. It is also interesting 
to note that the qualitative features of the phase dia- 
gram are very sensitive to changes in the next-nearest- 
neighbor (NNN) hopping constant t', in qualitative agree- 
ment with the unrestricted HF results. 12 In particular, 
while a pure AF phase is possible only at half filling for 
t' = 0, the AF phase can be doped by electrons, but 
not holes, for t' < at larger values of U, in agreement 
with previous results obtained with a more complicated 
method. 13 

The plan of the rest of this paper is as follows. In 
the next section we explain and justify our method using 
physical arguments. Mathematical details can be found 
in Sec. III. Section IV contains our conclusion and a sum- 
mary. 



II. THE METHOD 

We now explain our method, concentrating on the 
point where we deviate from standard MF theory. Precise 
mathematical formulas implementing this method will be 
given in the next section. As a representative example 
we discuss the computation of the phases by our method 
for the 2D Hubbard model with U = 6, and t = 1 and 
t' = —0.16 (see Eq. © below for the precise definitions). 
One reason for this choice is that it shows nicely several 
qualitative features which can occur in the phase dia- 
gram, another that these parameter values are of interest 
for high-T c compounds. 14 

MF theory for the Hubbard model is obtained by re- 



FIG. 2: Phase diagram of the 2D Hubbard model as a function 
of U and doping p for parameters t — 1, t' = —0.35, L = 60 
and j3 = 1000, computed as Fig. 1. For large U and p close to 
zero it becomes numerically difficult to distinguish between 
the F and P phase, which is the reason for the fuzzy phase 
boundaries in this region of the phase diagram. 



stricting HF theory to translational invariant states de- 
scribing antiferromagnetism (AF), ferromagnetism (F) 
and paramagnetism (P)££ It would be straightforward 
to generalize this and also allow for charge-density waves, 
ferrimagnetism etc. One thus starts with three varia- 
tional states which all are Slater determinants 1 ^ built of 
one-particle wave functions which are eigenstates of a 
mean field Hamiltonian where the Hubbard interaction 
is replaced by external field terms, 



I Slater) = IX), X = AF, F or P. 



(1) 



These fields include the the fermions density p and the 
magnetization which is staggered for AF, constant for F, 
and zero for P, and they are determined by the usual 
Hartree-Fock equations. It is important to note that the 
fermion density is fixed in the standard Slater states, but 
we use a generalization of Slater's variational principle to 
Gibbs states allowing for finite temperature and where 
the fermion density is varied by changing a chemical po- 
tential /i (grand canonical ensemble). 6 ' 16 "^ We now com- 
pute the Hartree-Fock ground state free energy per site, 
J~x, for each of these states X =AF, F and P, as a func- 
tion of /J,. 

Figure 3 gives the result for our example. At fixed 
value of /j,, the mean field ground state is determined by 
the minimum, 



T u 



X=AF,F,P 



Tx- 



(2) 



It is now important to note that the fermion density can 
be computed as derivative of the free energy as follows, 
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dT n 
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FIG. 3: Mean field free energy Tx of the 2D Hubbard model 
with t = 1, t' = -0.16, (7 = 6, L = 60 and /3 = 1000 as a 
function of the chemical potential p. Shown are the curves for 
X =AF, F and P (thin lines) and the absolute minimum !F m i n 
(thick line). The dashed lines indicate the particular values 
fii, i — 1,2, 3, of p where the phases change. At these values 
the derivative of .Fmm has discontinuities, and this leads to 
doping regimes with mixed phases; see Fig. 4. 



we use conventions such that particle-hole symmetry is 
manifest, p — 1 —* 1 — p corresponds to p — > —p and 
t' — ► — £'. From Figs. 3 and 4 is is obvious that this 
function p — 1 is, in general, only piecewise continuous, 
and it has jumps at the particular values of p where the 
minimum free energy curve changes, for example, from 
the AF to the F curve at the value p = pi. The physical 
interpretation of this is as follows. We start at p = 
where we obviously have the AF ground state and half- 
filling, p — l = 0. As we decrease p, p — 1 remains zero 
since .Faf does not change. This is due to the AF gap: as 
long as p remains in the gap the fermion density cannot 
change. For large enough p values the AF band edge is 
reached and the slope of J~af starts to decrease. However, 
before this can happen the F free energy has become 
lower and taken over: as one decreases p the F free energy 
decreases, and at a value p = p2 the two curves cross, 
Taf = Ff at p — pi- At this point we go from the AF 
to the F phase. Since the fermion densities px{pi) — 
1 = —dJ-x/dp\^ =il2 for the states X =AF and X =F 
are different, it is impossible to get a density value in 
between with either state. There is, however, a possibility 
to realize such a fermion density with the following state 
exactly at p = pi, 



|mixed) = w\AF) + (1 - w)|F), 



(4) 



with the relative weight w determined by the density as 
follows, 

p = wpAv{p2) + (1 - w)p F (p 2 ), 0<w<l. (5) 



FIG. 4: Doping p — 1 of the 2D Hubbard model as a function 
of the chemical potential p. The parameters are as in Fig 3 
(i = 1, t' = -0.16, U = 6, L = 60 and /3 = 1000). The 
curves are the derivatives of the corresponding ones in Fig. 3. 
The thick line determines the mean-field phase diagram, with 
the discontinuities at p — Pi, i — 1,2,3 determining doping 
regions where no pure phase F, AF or P is thermodynamically 
stable. The wiggles of the curves are due to finite size effects 
which, however, have no effect on the phase boundaries (this 
is demonstrated in the inset of Fig. 5). 



We now discuss the interpretation of this mixed solu- 
tion. One possibility is that the system has phase sepa- 
rated and split up into AF and F regions. 1 - Of course, the 
spatial structure of the actual state is not available in the 
MF description by the mixed state, but it can in principle 
be calculated using unrestricted HF. However, since the 
bulk free energy dominates over the interfacial free ener- 
gies in the thermodynamic limit, the mixed state gives an 
accurate description of the thermodynamics. We stress 
that the appearance of such a mixed state does not nec- 
essarily mean phase separation. The effect of the phase 
boundaries and other possible states have been excluded 
in our approximation. To know the actual state in the 
mixed regions thus is beyond our calculation and can be 
decided only by doing more work, e.g., using unrestricted 
HF taking into account more complicated states. Nev- 
ertheless, the occurrence of such a mixed states proves 
that no simple translational invariant state of the kind 
assumed in our MF ansatz is thermodynamically stable. 
The mixed regions of the phase diagram are of particu- 
lar interest since there the free energy is degenerate and 
thus the details of the solution can be strongly affected 
by fluctuations, phase boundaries, or details neglected in 
the model. 

It is important to note that there are two further jumps 
of p and two further corresponding mixed phases: one at 
p = pi with F coexisting with P, and another at p = p% 
with AF and P coexisting. It is also interesting to note 



that, while for t' = the mean field free energies are in- 
variant under the electron-hole transformation ji — » —fi, 
the finite value of t' = —0.16 here leads to a qualitative 
difference between hole doping (/i < 0) and electron dop- 
ing (/i > 0). As seen in Fig. 3, the F state can compete 
with the AF state only for \i < 0, and this implies that 
it is possible to dope the AF state by electrons but not 
by holes. 

We thus see that, even though we restricted Hartree- 
Fock theory to simple translation invariant states as in 
Eq. ifTjl, our way of treating the doping constraint has 
implicitly also included the possibility of having a mixed 
state as in Eq. J3J as groundstate, and we find that such 
a mixed state indeed occurs in a significant part of the 
doping regime. 

We stress that our method to determine the phase 
boundary does not increase the computational effort of 
mean field theory, and it is easy to do the computations 
also for large system sizes. Most of our computations 
were done for a L x L lattice with L — 60. While at this 
values of L some finite size effects are still visible in the 
relation between doping p and the chemical potential /i 
(see Fig. 4) , the inset in Fig. 5 demonstrates that result- 
ing phase boundaries are practically identical with the 
ones in the thermodynamic limit. We also checked that 
the value [3 — 1000 we used for the inverse temperature 
practically gives the zero temperature phase boundaries. 



III. FORMALISM 

We now give the formal implementation of our method. 
We start by fixing our notation. We consider the 2D 
Hubbard model defined by the Hamiltonian 

H = -* E c U c i,a - *' E c U c j,« + H - c - 



{i,j),a 



((*J»>« 



'-E c L c v* + C/ E^T ~ o)Kl _ o) 



(6) 



with the on-site repulsion U > and the hopping ampli- 
tudes t > and t' between the nearest neighbor sites (i, j) 
and next-nearest neighbor (NNN) sites ((i,j)) on a square 
lattice with L 2 sites, respectively; the fermion operators 
c\ a are parameterized by the spin variable a =\, | and 



lattice sites i 



(i x ,i y ) where i x>v 



1,2, 



, L, and 



c] a c i a are number operators, as usual. The 



fermion density is 



^Ev^.«> 



V- 



(J) 



with (•) the ground state expectation value to be specified 
below. 

We recall that unrestricted Hartree-Fock (HF) theory 
is formally obtained by introducing 



q % = (m), m, ; = (s^ 



(8) 



=> 10 




FIG. 5: Phases of the 2D Hubbard model as a function of the 
chemical potential for the same parameters as in Fig. 2 (t = 
1, t' = -0.35, /? = 1000, L = 60). Inset: Blowup of the region 
around the minimum of the phase lines in the main figure, 
showing interesting fine structure in the phase diagram. Also 
shown is the result from a calculation for system size L — 120 
(crosses). The coincidence between results for two different 
system sizes demonstrate that L = 60 is practically already 
in the thermodynamic limit. 



and replacing the Hubbard interaction by external field 
terms as follows, 



1 



n i,T n U 



K 



pf) + \{^ 



where ^[/m^ and -^Uqi are mean fields coupling to the 
fermion spin Sj = Y^ a ,a' c i,a flr aa' c i,a' and (local) den- 
sity rii = J2 a c l, a c i, a i respectively; a = (ci, ct 2 , cr 3 ) are 
the usual Pauli spin matrices. This replacement leads 
to a Hamiltonian describing non-interacting fcrmions in 
external fields, H — > Hhf with 



U 
1 



H HF = J2^(™ 2 i-qf)+ E C Xcc\ a ,i, a -% a ' ( 9 ) 



where 



/^i,a, 



-UjSaa 1 + Sijl—U[. 

(qi - l)S aa 



m,; • (T aa i + 



(10) 



is a self-adjoint 2L 2 x 2L 2 -matrix which can be inter- 
preted as a one-particle Hamiltonian. One now interprets 
(•) in Eq. JSJ) as the expectation value in the ground state 
of Hhf in Eqs. [J9I10J1 . This yields the HF equations al- 
lowing to self-consistently compute qi and m^ (see e.g. 
Sec. II in Ref. El. 



We now observe that these HF equations can also be 
obtained as saddle point equations dT jdva.i = dT/dpi = 
from the free energy function 



T - 



1 



2 \ogZ 



where 



Z = Tr 



-PH h 



(11) 



(12) 



is the partition function defined by a trace over the 
fermion Hilbert space, and (3 is the inverse temperature. 
A straightforward computation yields 



T 



U 
AL 2 



E(* 



2L- 



P?) 



1 \S v,^ 

2 log cosh — — . 



L 2 



(13) 



with Eg the eigenvalues of the one-particle Hamiltonian 
h = (hi >a ,j,a') in Eq. [fTUjI . 6 

The physical solution of the HF equations are such that 



T a 



min max .F(mj, <# 

mi qi 



(14) 



see Ref. [IT] for a mathematical proof or Ref. |6| for a 

derivation using functional integrals. The correspond- 
ing fermion density is then given by Eq. 10. We stress 
that Eq. I|14|) , while implying standard HF theory, is not 
equivalent to it: the standard HF equations can have sev- 
eral solutions, but Eq. (|14|) provides a simple method to 
solve HF equations so as to avoid the unphysical solu- 
tions: first maximize T with respect to the qi, and then 
minimize with respect to the m^. In case we restrict HF 
theory by making a simplifying ansatz for the mean fields 
qi and m.i as below, it can happen that one finds several 
HF solutions at a fixed value of p. In this case one must 
take the solution minimizing T . 

Mean field theory is obtained from HF by restricting 
to mean fields which are invariant under translations by 
two sites. For the different states discussed in this paper 
one further simplifies to 



(15) 



where e z is the unit vector in ^-direction. With this 
restrictions it is easy to compute the eigenvalues Eg by 
Fourier transform. One obtains 



AF 


q t = q, 


m, = m AF (-l) lx+iH e 


F 


it = q, 


u\i = mFe z 
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Qi = q, 


mj = 



AF: 



F 
P 



Ei 



k.± 



i[e(k)+e(k + Q)] 



+ U{q - 1) - M ± 



±^/[e(k)-e(k + Q)]+([/m A F) 5 
E K ± = e(k) +U(q-1± m F ) - p, 
E K ± = e(k) + U(q-l)-n 



(16) 



where the quantum numbers labeling the eigenvalues are 
£ = (k, e) with e — ± a band index and k = (k x ,k y ) 



with k x , y = (2tt/L)x integer momenta restricted to the 
Brillouin zone — 7r < k x ^ y < n; Q = (jr, w) is the AF 
vector, and 

e(k) = — 2t[cos(k x ) + cos(fc y )] — At' cos(k x ) cos(k y ) (17) 

is the usual tight binding band relation. Thus the mean 
field free energy becomes 



U 

1 



Fx = -r( m x -q 2 ) 



"72 E co^o^k 



k,e=± 



(18) 



for X =AF, F and P (rnp = 0), where the k-sum becomes 
an integral in the thermodynamic limit L — ► oo. The 
standard mean field equations (see e.g. Sec. II in Ref. 4) 
are obtained from this from differentiation, dTx/dq = 
dFx/dmx — 0. Note that q = px (fermion density at 
fixed /i in the A-state) but, as explained in Sec. II, the 
relation of px to the system density p is somewhat subtle. 



IV. CONCLUSIONS AND DISCUSSION 

In conclusion, we have presented a simple generaliza- 
tion of standard mean field theory, including the possi- 
bility of phase separated mean field states. We have pre- 
sented results for the phase diagram of the 2D t — t' — U 
Hubbard model, including values of parameters sug- 
gested by the high-T c materials. We find that the NNN 
hopping t' significantly alters the solution. The result- 
ing rich and nontrivial phase diagrams show significant 
qualitative differences between electron and hole doping. 
Moreover, a finite t' suppresses order in the weak cou- 
pling regime, but can have the opposite effect at strong 
coupling; see Figs. 2 and 5. Thus the results presented 
here are much richer than those obtained by standard 
MF theoryii4 The correctness of our method is justified 
by mathematical rigorous results^ 

We stress that the method presented here does not 
necessarily produce accurate solutions to the problem, as 
is often the case with mean field theory. Nevertheless the 
method provides a useful starting point for estimating the 
structure of the phase diagram, providing cheap guidance 
for more accurate but costly calculation methods towards 
interesting regimes in the phase diagram. 

The simple theory presented here can be straight- 
forwardly generalized to a number of interesting cases, 
including more general mean field states like ferrimag- 
netism or stripes, and to more complicated models with 
additional interaction terms or more bands, etc. 
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